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The dynamics of a galactic disk in a non-axisymmetric (triaxial) dark halo is stud- 
ied in detail using high-resolution, numerical, hydrodynamical models. A long-lived, 
O ■ 

' two-armed spiral pattern is generated for a wide range of parameters. The spiral 

> ■ 

O ' structure is global, and the number of turns can be two or three, depending on the 

g _ 

model parameters. The morphology and kinematics of the spiral pattern are studied 

OO ■ 

I as functions of the halo and disk parameters. The spiral structure rotates slowly, 

' and its angular velocity varies quasi-periodically. Models with differing relative halo 

o' 

^ . 

Q_|, gaseous component, and angular rotational velocities of their halos are considered. 

o 

^ ■ 1. INTRODUCTION 



masses, halo semi-axis ratios, distributions of matter in the disk, Mach numbers in the 



Galactic gaseous and stellar disks are axisymmetric only in a first approximation. Spi- 
ral density waves, a central bar, the non-axisymmetry (triaxial nature) of a dark massive 
halo, satellite galaxies, or a triaxial bulge can disturb the axial symmetry of the "external" 
gravitational potential in which the galactic disk component is situated. Results of dynamic 
modeling of specific disk galaxies, designed to agree with photometric and kinematic data, 
indicate the need to take into account a dark halo of mass Mh (which exceeds the mass of the 
^ ' disk component Md by a factor of one to four [1-5]) inside the optical radius of the stellar 

component R. The kinematics of objects in the immediate vicinity of the MilkyWay-M31 
pair indicate that the total mass inside 0.96 Mpc is ~ 2 x IO^^Mq [6], and the ratio 
of the masses of the MilkyWay and M31 is Mmw/Mmsi = 4/5. Since the mass of the disk 
component of our Galaxy probably does not exceed 5 x 10^^ Mq [7], ~ 15% of the total mass 
within 1 Mpc belongs to the disk subsystems. 

There is much theoretical and observational evidence indicating the absence of central 
symmetry in the dark-mass distribution within the disk components of S galaxies. The 
triaxial shape of the halo can be characterized by three space scales: a, b, c. 

Studies of galaxy formation based on cosmological N-body models indicate triaxial mass 
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distributions in cold-dark-matter (CDM) halos, with the following semi-axis ratios obtained 
in different studies: s ^ c/a ^ 0.6 - 0.7, g = 6/a ~ 0.8 - 0.95 [8]; (s) ~ 0.7, (g) ~ 0.8 
[9]; (g) ~ 0.75 [10]. At the same time, triaxial parameters s and q depend on the halo 
mass, redshift, and distance to the halo center. The dark halo is more symmetrical near 
the center than at the periphery. The distributions turn out more symmetrical in hot-dark- 
mattermodels than in CDM models [9]. 

Observations of the Sagittarius tidal stream (Sagittarius dSph) indicate a lack of cen- 
tral symmetry of the halo, with the axial ratio c/a >0.8 [11]. However, if this stream is 
dynamically young, then cja c^. 0.6 is possible [12] . Arc-shaped stellar tidal streams found 
in the Sloan Digital Sky Survey (SDSS), apparently resulting from the absorption of dwarf 
galaxies, seem to be promising indicators of a dark halo. The Virgo over-density structure 
could be related to such peculiarities in the distribution of the observed matter in the halo 
[13]. Stellar tidal streams have been discovered in M31 [14] and a number of outer galaxies 
(e.g., NGC 5907 [15], NGC 3310 [16], NGC 1055 (see http://www.imagingdeepsky.com)). 
These streams are apparently connected with the disruption of dwarf galaxies in the massive 
triaxial halo. 

Observations of stars (RR Lyr, F dwarfs) indicate a lack of central symmetry of the 
halo [13]. Direct kinematic observations of various halo objects (globular clusters, satellites, 
blue horizontal-branch stars, red giants) can be used to find the radial velocitydispersion 
profile which, in turn, indicates the character of the gravitatingmass distribution for a given 
halo model. New possibilities for deriving halo parameters are provided by observations of 
hyper-velocity stars [17-19]. 

Deep surveys of edge-on galaxies indicate a strong prevalence for stellar halos in the form 
of oblate spheroids with axial ratios of ~ 0.6, which have been attributed to tidal interactions 
[14, 20]. Cosmological models taking into account the stellar-halo dynamics also give triaxial 
shapes for the dark halo (s = 0.84, q = 0.94) and stellar halo (s = 0.84, q = 0.91) [21]. 

Several techniques for determining halo shapes (a, ft, c) from observations can be identified 
[22], based on: (1) expansion of the HI gaseous disk with distance from the center (flaring of 
the gaseous layer) — other conditions being equal, the thickness of the gaseous disk decreases 
with increase in the degree of flattening of the halo; (2) warping of the gaseous layer; (3) 
X-ray isophotes; (4) outer polar rings— a dynamical model of the polar ring is constructed 
allowing for an asymmetric halo and fit to the observational data; and (5) the parameters 
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of a precessing dust disk. The observational data for our Galaxy indicate that the both the 
stellar and dark-matter halos do not display central symmetry. Different approaches based 
on observations of different objects lead to estimates of s ~ 0.2 — 0.9 (see, e.g., [23, 24]). 

Thus, one of the most important properties of a halo is the triaxial character of the mass 
distribution in it. In turn, information on the character of the darkmatter distributions in 
galaxies and in their immediate vicinities, together with observations of the isotropic gamma- 
ray background, can be used to derive constraints on particle dark-matter models [25]. 

Let us consider the dynamics of a rotating gaseous disk in the non-axisymmetric gravita- 
tional field of a dark halo. Let us assume that themass of the stellar-gaseous disk within the 
limits of the optical radius R is comparable to the dark mass [2, 26, 27]. The influence of 
the global asymmetry of the potential on the formation of dynamical structures in the disk 
is of the most interest. It seems natural to expect that the disk responds to this influence 
by forming a spiral pattern [28, 29]. 

Let us consider the hypothesis that the nonaxisymmetry (triaxial nature) of the dark halo 
is responsible for the formation of the spiral structures in galactic disks, at least for some 
galaxies. We will restrict our consideration here to the gaseous disk. A separate study will 
be devoted to taking into account the stellar disk. 



2. MODEL OF A GASEOUS DISK IN THE EXTERNAL POTENTIAL 

The character of the spatial mass distribution in the dark halo determined from observa- 
tions is uncertain. Several models for the halo density Qh have been proposed, such as the 
models of Navarro et al. [30] (NFW) and Burkert [31], the exponential model [32], and the 
quasi-isothermal halo model [1, 33] (iso). We will generalize this last model for the case of 
a triaxial halo, by writing the following expression for the potential 

vl/,(x, z) = AnGg,oa' ■ |ln(C) + + ^ In | , (1) 

— H — r- H — r-, and a^, a„, and are characteristic scales along the correspond- 
ed fly ai 

ing axes. The Poisson equation gives the density distribution 

arctan(^) 
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In the limiting case of a centrally symmetric halo a — ax — ay — az, have a quasi- 
isothermalmodel halo with the density distribution 

Q{r) = . j!'^ ^2 > (3) 

1 + [r/ay 

from (1), (2); this density distribution ensures constancy of the rotational velocity at large 
distances r ^ a, where the rotation curve has a plateau. The halo models differ in the circular 
velocity Vc at large distances (r >-Rmax, -Rmax is the boundary of the stellar disk) [34]. In 
the inner region of the galaxy, we have Vj''°^ < ^/"^ < Vf^^ [34]. Therefore, other things 
being equal, our choice of the model (1) should be considered to yield a lower limit for the 
influence of the non-axisymmetric halo. 

Below, we restrict our consideration to the model with ^ ay = a^. If the halo is 
flattened in the vertical direction (a^ < a^^y), this effectively increases the halo mass, since 
the gravitational force of the halo in the z — plane becomes greater, so that the case 
a^ — Ox again gives us a lower limit. When q — a^/ay 7^ 1 in the plane of the disk, we 
have a non-axisymmetric halo that can, in general, rotate with a small angular velocity Q^. 
Observational data and cosmological models allow values eh = lo' — 1| ~ — 0.2 [8-10]. 

The character of the flow in the gaseous disk depends considerably on the values of the 
indicated model parameters. However, the generation of spiral structure is observed in all 
cases, and we have a longlived, quasi-stationary, two-armed wave in the main region of the 
disk (0.5 r < 2). Figure 2 shows the typical long-term evolution of the disk [t — 40 
corresponds to six rotations of the disk periphery r ~ 2). The two-armed perturbation 
covers almost the entire disk, and is trailing except at the center itself, where more complex 
structures are observed. 

Let us assume that the rotational angular velocity of the halo 0/j is less than the rotational 
angular velocity of the outer edge of the disk, VL[R) — V{R)/R. If typical values of the linear 
velocity V{R) for galactic disks in the plateau region are Vmax — 100 — 300 km/s at distances 
i? ~ 10 — 15 kpc {Q{R) ~ 10 — 30 km s~^kpc~^), this imphes for the angular velocity of the 
halo Qfe <5 km s~^kpc~^. 

Let us consider the complete system of gasdynamical equations in the form 

^ + V.(pu) = 0, (4) 
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dpu 



+ V • (puu^ + PI) = -pV^ 



(5) 



dt 
dE 
'dt 



+ V • ([£; + p]u) = -pvL ■ V* , 



(6) 



where t is the time, p the gas density, P the pressure, u = [u^v^wY the velocity vector, 
/ = diag[l, 1, 1] a 3 X 3 unit tensor, and = + the gravitational potential, composed 
of an external component (halo + stellar disk) '^q{y, t) and the potential of the selfgravitating 
gas ^g{r, t). The total energy per unit volume of the gas is 



where e is the specific internal energy. The equation of state of the gas is defined by 
e = P/(p(7 — 1)), where 7 is the adiabatic index. The system of equations (4)-(7) must 
be supplemented with the Poisson equation = AnGg, to take into account the self- 

gravity. 

Wc used total variation diminishing (TVD) version of the MUSCL scheme [35, 36] to 
numerically model the system of gas-dynamical equations (4)-(6). We approximated the 
equations using the finite-volume method, which ensures that conservation laws are obeyed 
at the discrete level [37]. This numerical scheme is related to Godunov schemes with sec- 
ondorder accuracy in time and third-order in space. The modified HLLC method [38] was 
used to calculate the fiuxes of physical quantities through cell boundaries; this provides a 
means to model shocks and contact and tangential discontinuities through the boundaries. 
The construction of the scheme with second-order accuracy in time was carried out using 
a predictor-corrector scheme compatible with the condition that the total variation of the 
numerical solution not decrease (the TVD condition [36]). A third-order approximation in 
space was used in regions of smooth fiow, and was achieved using the MUSCL procedure for 
reconstructing simple quantities (the density, pressure, and velocity vector) at the bound- 
aries of the computation cells. This simple and computationally efficient method makes it 
possible to obtain shocks spread over three cells. The corresponding parameter responsible 
for the artificial compression was used to improve the resolution of tangential and contact 
discontinuities in the interpolation. 

The computations were carried out on cylindrical (r, </?, z) and Cartesian {x, y, z) grids 
with high resolution. The dimension of the numerical grids for the two-dimensional (2D) 
models reached 5000 x 5000 for the Cartesian grid and {N^ x N^) = (2400 x 720) (for the 




(7) 
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polar grid. The grid parameters in the three-dimensional (3D) model were Nr — 600, — 
360, Nz — 200. A special numerical algorithm was constructed for the through calculation of 
the gas-vacuum interface, which preserved the positivity of the scheme [39] . In this case, the 
matter is situated within the computation domain. A comparative analysis demonstrated 
the extremely small influence of the boundary conditions on the evolution of the constructed 
models for the gaseous disks. 

In order to estimate the influence of the numerical technique used on the disk dynamics, 
we also carried out hydrodynamical calculations using Smoothed Particle Hydrodynamics 
(the SPH method) [40, 41]. The SPH algorithm is related to fully Lagrangian algorithms, and 
so is free of the need to use the approximation of a spatial grid; this removes a considerable 
number of theoretical and algorithmic difficulties. The scheme is constructed based on the 
motion of particles whose evolution in time and space directly reproduces the conservation 
of mass, momentum, and angular momentum. A modification of the standard predictor- 
corrector scheme was used for the time integration of the SPH equations; this modification 
ensures second-order accuracy in time. The calculations were carried out N — 10^ — 10^ 
particles, and demonstrated good agreement with solutions based on the grid TVD method. 
The selfgravity was taken into account using the TREEcode algorithm. 

The initial equilibrium state of the gaseous disk is defined by the radial and vertical 
profiles of the rotational velocity V, density g, and sound speed Cg in the absence of radial 
and vertical motions, u — w — 0. It is convenient to characterize the state of the disk using 
the Mach number, M. — V/cg. Since M. depends on the coordinates, it is convenient to 
introduce the quantity Aio = ymax/cs(r = 0). Figure 1 presents typical radial profiles of the 
disk parameters at the initial time for the assumed dimensionless quantities. The rotation 
curves are typical for galactic-disk components, and are characterized by a small section with 
a quasi-rigid-body rise at the disk center and a more extended section with roughly constant 
velocity (the so-called plateau in the rotation curve). 

The numerical model assumes G — 1, and the galaxy mass within the hmits of the radius 
r < 1 is Mgai = 4. The optical radius of the galaxy (the outer edge of the stellar disk) is 
i? = 2 in dimensionless units. In this case, a typical value of the dimensionless rotational 
velocity of the disk is \/jnax — 2, and the rotational period of the periphery of the gaseous 
disk (r ~ 2 — 3) is T ~ 6 — 10. The gaseous disk extends beyond R. 

We considered various initial radial profiles of the gas surface density a: 
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(1) (T oc exp(— r/L), 

(2) ao,{l-r{rlUf-r{rlL^f)-\ 

(3) a«(l + (r/L3n-5/2, 

(4) a oc cosh~^(r/L4). 

The calculation domains in the various models are situated within -Rmax = 5 — 15, where 
c(-Rmax)/c(0) = 10~^ — 10~^°, which ensures an absence of influence from boundaries. 

At the initial time t = 0, the axisymmetric equilibrium disk is in the axisymmetric 
potential ^{r,t — 0) with Sh — ^ {o-x — o-y)- Then, the value of the scale was increased 
linearly to — 0.001 -j- 0.2 over a time Th for the various models. After the nonaxisymmetry 
was "switched on" at times t > Th, the halo potential remained constant. 

In the following section, we describe the main results of our modeling of the disks, fo- 
cusing on the possible generation of various structures by the triaxial halo, considering the 
morphology and kinematics of these structures based on the results of over 80 numerical 
experiments. 

3. GASEOUS DISK IN A NON-AXISYMMETRIC HALO 

The main free parameters of the constructed model and their typical dimensionless values 

are: 

(1) the degree of non-axisymmetry of the halo 6h = 0.2 (g = ay/a^ = 1 — Sh)- 

(2) the width of the quasi-rigid-body section of the rotational velocity in the central region 
of the disk d = 0.015 - 0.5; 

(3) typical scales for the radial distribution of the surface density L = 0.5 — 1.5; 

(4) the sound speed of the gas or the disk temperature, characterized by the effective Mach 

number A^q = 5 — 30; 

(5) the relative gas mass Hg = Mg/M^ (Mq is the galaxy mass within the radius r < R = 2); 

(6) the rotational angular velocity of the halo flh — — 0.2, corresponding to the corotation 
radius at the periphery of the gaseous disk; 

(7) the transition time from the axisymmetric to the non-axisymmetric halo — — 30. 
Let us note some typical features of the morphology and kinematics of the disk indicated 

by our numerical simulations: 
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1. In all cases, the presence of a non-axisymmetric halo leads to the formation of spiral 
structures in the gaseous disk (see the distribution of the logarithm of the surface 
density in Figs. 3a-3d). Under typical conditions, perturbations grow to the strongly 
nonlinear stage over one to three disk revolutions, forming a system of shocks (Fig. 3e). 
The time when the quasi-stationary spiral structure forms increases with decreasing 
Eh- With decreasing the spirals become more tightly wound and have smaller 
amplitudes, other conditions being equal (Fig. 3e). A massive, non-axisymmetric 
halo can generate nonlinear waves even when Sh < 0.01, but the rise time in this 
case exceeds eight rotation periods of the disk periphery. The formation of the spiral 
structure in the essentially nonlinear stage depends only weakly on the time when 
the non-axisymmetric part of the potential is switched on, t/j, which was varied from 
Th — (instantaneous) to an adiabatically slow ramping up during several rotations of 
the disk periphery (about 1-1.5 billion years for a typical S galaxy). 

2. The analysis of the gas-flow structure indicates a complex character that includes a 
system of shocks and of regions of strong shear flows (Fig. 4). The nonlinear (shock) 
waves have large-scale structure, which encompasses the disk right to small densities 
q/ Qmax ~ 10~^. The shape of the spiral depends substantially on the model parameters, 
but it is a large-scale formation, formed by rotating through the angle of ~ 27r, or, in 
some cases, by more than An. The amplitude of the spiral waves depends on the degree 
of non-axisymmetry of the halo Eh, but the formation of the shocks occurs for all of 
the parameter values considered. A zone of shear flow is formed in the region of the 
spiral density wave near the shock front (Fig. 4a). The region of strong variation of 
the tangential velocity component is large compared to the shock front, whose width 
i is determined by the properties of the numerical model (Fig. 4b). The quantities a 
and u exhibit a jump at the shock front (within i), in contrast to v, whose region of 
variation is substantially larger than £. The width of the shock front £ decreases with 
increasing resolution of the numerical scheme (increasing numbers of cells Nj., N^), so 
that the shock front is described by three to four numerical cells; this was achieved 
by using a numerical scheme with second-order of accuracy in time and third-order 
accuracy in the spatial coordinates. The width of the shear flow region is flnite and 
does not depend on the choice of parameters for the numerical grid. 
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3. The geometry of the spiral patterns depends on the gas- rotation curve V{r), the char- 
acter of the density distribution in the dark halo gh{r,(p), and the radial profiles of 
the disk surface density (T(r)and the sound speed Cs{r). Leading spirals, which always 
become trailing spirals at the periphery and produce complex structures (embedded 
structures are possible) at the center, can be formed in the central part of the disk 
(Fig. 5) . In all the figures presented here, the gaseous disk rotates in the counterclock- 
wise direction. Thus, trailing spirals untwist in the clockwise direction, and leading 
spirals in the counterclockwise direction. The presence of leading spirals in the very 
central zone of the disk is clearly visible in Figs. 3d, 5c and Figs. 9a, 9b, 10a. 

4. The generation of pronounced structures formed by leading spirals in the central 
region occurs in gaseous disk models without allowance for the perturbation of the 
stellar-subsystem?s density. Since structures are not formed in models with coUi- 
sionless stellar disks, the central structures will be weaker in self-consistent stellar- 
gaseous systems, and the question of the dynamics of the central region requires a 
separate study. Note that complex gaseous structures manifest as, e.g., dust streaks 
inside the bar are clearly visible in images of SB galaxies (e.g., NGC 1672, NGC 1300, 
NGC 1365, NGC 3627, etc.). In addition to gaseous?dusty spirals emerging from the 
center, on the whole, along the major axis, structures tracing the shape of the bar 
itself (e.g., in galaxies NGC 1097, NGC 1288, NGC 2903) or outer spiral arms that 
do not start from the bar ends but encompass part of the leading edge of the bar 
(NGC 4548, VCC 1615) are observed. The socalled three-kpc arms in our Galaxy may 
serve as an example. 

5. As in the case of decreasing Sh, decreasing the zone of the quasi-rigid-body rotation in 
the center of the gaseous disk promotes the formation of more tightly wound spirals. 
Simultaneously, the conditions for the generation of structures in the center become 
less favorable. Trailing spirals now arise in the center itself and encompass the entire 
disk (Fig. 5). 

6. In cooler disks (with higher Mach numbers), more tightly wound and thinner spirals 
are formed in the nonhnear stage (Fig. 6). 

7. The typical spiral pitch angles are i ~ 5° 20°. More open patterns do not formin 
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themodels considered. Figure 7 shows an example of an approximation of the obtained 
structure with sections of logarithmic spirals r = raexp{x^) with pitch anglesi = 
7r/2 — arccos(x/A/l + X^)- For the case presented, i increases from 10° in the central 
zone to 14° at the periphery. This small increase in i with radius is typical for the 
models constructed. 

8. Let us now consider the rotational angular velocity of the spiral pattern Qp. Figure 
8 shows the time dependences of Qp for the two-armed spiral at various radii r. An 
initial stage in the pattern formation can be distinguished (t < 10; these processes 
proceed faster nearer the center than at the periphery), after which a quasi-periodic 
rotation regime for the non-linear perturbations develops. We emphasize that the spiral 
pattern rotates even in the case of a stationary halo — 0. However, the rotational 
angular velocity is then low, and the corresponding corotation radius is beyond the 
disk. Computations with a slowly rotating halo, Qh ^0.2, show that the geometry and 
kinematics of the spiral patterns depend weakly onQh. A small increase of the angular 
velocity of the wave in the disk is observed with increasing fl^ ^. Objects whose spiral 
structure rotates fairly slowly and whose corotation radius is at the periphery of, or 
even outside, the outer hmits of the spiral structure (e.g., in NGC 3359 [42]) are of 
interest for the spiralstructure formation mechanism discussed here. 

9. There is qualitative agreement between the results of the 2D and 3D gas-dynamical 
computations (Figs. 9, 10). In the initial stages of formation of the spiral pattern 
(over several rotations), the formation of strongly nonlinear wave fronts propagating 
in the vertical direction at an angle to the plane of the gaseous disk is observed. The 
geometry of the trailing spiral pattern in the main part of the disk differs slightly in 
the 2D and 3D models (Figs. 9, 10). Discrepancies are observed only in the central 
region, where the formation of complex, small-scale structures is possible. In this case, 
the presence of a vertical dimension can substantially change the gas dynamics, since 
the sound speed in the center is known not to be small compared to the rotational 
velocity of the system, and vertical motions are of considerable importance [43]. 



Note that the question of the location of the corotation radius is far from resolved, and this is one of the 
key questions in spiral structure theory. 
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10. On the whole, taking self-gravity into account increases the efficiency of the formation 
of the spiral pattern by the non-axisymmetric halo. We have considered only grav- 
itationally stable disks. A change of the wave geometry is observed, and the wave 
amplitude increases (Fig. 11). 

In conclusion, let us summarize the main differences of the spiral-pattern modeling in a 
nonaxisymmetric halo and in the presence of a central bar: 

(1) the bar is a rapidly rotating structure relative to the spiral pattern; 

(2) the dimensions of the central stellar bar are small compared to those of the galactic 
disk, and the effect of the bar is dominant at the center and rapidly decreases with radius; 
the triaxial halo influences the entire disk, and influences the disk center only to a smaller 
degree; 

(3) the degree of non-axisymmetry of the halo is substantially less than of that the bar, 
so that the geometry of the spiral waves differs considerably from the disk in the nonax- 
isymmetric halo, even inside the bar, where the gas, on average, rotates faster than the 
bar. 

The results of our numerical modeling show that a triaxial (non-axisymmetric in the plane 
of the disk) halo can lead to the generation of a spiral pattern in the galactic disk, whose 
properties depend strongly on the parameters of the disk and halo. 

The question of whether non-axisymmetric halos are responsible for the real spiral pat- 
terns observed in galaxies is, of course, far from being resolved. The real picture can be 
made considerably more complex by the formation of transient spirals, a central stellar bar, 
or global modes [44] due to the development of gravitational instabihty. In any case, the 
models considered here can be used to place some constraints on the shapes of dark halos in 
the plane of their galactic disks. If the properties of the observed spiral pattern in a galaxy 
are incompatible with our modeling results, then this indicates a high degree of axisymmetry 
for the dark halo in the plane of the disk. 
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Figure 1. Typical dependences of tlie (1,2) rotational velocity V(r) (without and with a bulge, 
respectively), (3) surface density C7(r)/cr(0), and (4) sound speed Cs{r) on the radial coordinate for 

the initial equilibrium state of the gaseous disk. 
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Figure 2. Distributions of the logarithm of the surface density a in the plane of the disk at various 
times, for the model with d = 0.02, Sh = 0.1, Mq = 10, and L3 = 0.7. 




Figure 3. Distributions of the logaritlim of the surface density in the plane of the disk for models 
with Mo = 10,d = 0.1, and Sh = (a) 0.01, (b) 0.05, (c) 0.1, and (d) 0.15, all at the same time, (e) 
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Figure 4. Model with eh = 0.1, d = 0.02, L3 = 0.7, Mo = 10, and = 0. (a) Distribution of 
lg((divu)^) (in shades of gray), together with the velocity field u. (b) Radial profiles of the radial 
velocity u, azimuthal component v, and surface density a near the spiral density wave. The function 
values at the nodes of the numerical grid are marked near the curves. 




Figure 5. Spiral patterns with various geometries in the central region at the same time for 
A^o = 10 and £ = 0.1 , in the models with (a) d = 0.02, (b) 0.1, and (c) 0.3. 




Figure 6. Models with Macli numbers Mo= (a) 10, (b) 15, (c) 20, and (d) 30 witli £ = 0.1 and 

d = 0.1, all at the same time. 
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Figure 7. Geometry of spirals in the model with Ado = 10, £h = 0.1, and d = 0.02. Superimposed 
on the distribution of the logarithm of the surface density are curves showing sections of logarithmic 
spirals with various pitch angles: 1 i = 10°, 2 i = 12°, 3 i = 14°. 




Figure 8. Dependences flp{t) for the two-armed spiral for disk rings situated at different distances 

from the center r. 



Figure 9. Distributions of the surface density in the (a) 3D and (b) 2D models in the TVD 




Figure 11. Radial profiles of the surface density along the ray (p = in the SPH disk model (i) 

without and {2) with self-gravity. 



